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Abstract: Persian oak (Quercus brantii vat. persica) is a dominant tree species of Zagros forests in a semi-arid area, 
western Iran. However, the capacity of biomass and carbon stocks of these forests is not well studied. We selected 
three types of oak, i.e., seed-originated oak, coppice oak and mixed (seed-originated and coppice) oak of Zagros 
forests in Dalab valley, Ilam Province, Iran to survey the capacity of biomass and carbon stocks in 2018. Thirty plots 
with an area of 1000 m? were systematically and randomly assigned to each type of oak. Quantitative characteristics 
of trees, such as diameter at breast height (DBH), height, crown diameter and the number of sprouts in each plot 
were measured. Then, aboveground biomass (AGB), belowground biomass (BGB), aboveground carbon stock (AGCS) 
and belowground carbon stock (BGCS) of each tree in plots were calculated using allometric equations. The litterfall 
biomass (LFB) and litterfall carbon stock (LFCS) were measured in a quadrat with 1 mX1 m in each plot. One-way 
analysis of variance and Duncan's test were performed to detect the differences in biomass and carbon stocks among 
three types of oak. Results showed that AGB, BGB and BGCS were significantly different among three types of oak. 
The highest values of AGB, AGCS, BGB and BGCS in seed-originated oak were 76,043.25, 14,725.55, 36,737.79 and 
7362.77 kg/hm2, respectively. Also, the highest values of LFB and LFCS in seed-originated oak were 3298.33 and 
1520.48 kg/hm2, respectively, which were significantly higher than those of the other two types of oak. The results 
imply the significant role of seed-originated oak for the regeneration of Zagros forests. Further conservation strategy 
of seed-originated oak is an important step in the sustainable management of Zagros forests in Iran. 
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1 Introduction 


Increasing greenhouse gases including carbon dioxide (CO2) has been identified as the source of 
global warming or abrupt climate change (Allen et al., 2006; Lorenz and Lal, 2010). Human 
activities for reducing atmospheric CO2 were started when the impacts of global warming were 
specified (Sohrabi et al., 2016). The increase rate of atmospheric CO2 concentration can be 
reduced through the process of carbon sequestration (Lorenz and Lal, 2010). Carbon stocks in tree 
biomass and forest soils are economically considered as the most possible procedures to reduce 
atmospheric carbon (Brown, 1996). Trees play an important role in reducing atmospheric CO2 by 
absorbing and accumulating in their leaves, branches, stems and roots as biomass (McPherson and 
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Simpson, 1999). The aboveground biomass (AGB) in trees generally accounts for the greatest 
fraction of total living biomass in a forest (Brown, 2002). Forest productivity and aboveground 
carbon stock (AGCS) are also related with many factors, such as tree size inequality, stand age, 
nutrient regime and climate anomalies (Chen and Luo, 2015). Changes in stand characteristics 
through forest succession have significant impacts on forest productivity and AGCS (Becknell 
and Powers, 2014). This is because tree size inequality among and within species is critical 
toward maintaining the species diversity, diameter at breast height (DBH) and height (Wang et al., 
2011), which has been recognized to significantly affect forest carbon stocks (Lexergd and Eid, 
2006; Zhang and Chen, 2015). 

The amount of AGB in a region can be estimated either by direct or indirect methods. The 
direct method consists of cutting and weighing the AGB in an established area. This method is 
destructive and very time-consuming. Therefore, allometric equations are often used for 
estimating forest biomass. Allometric equations relate the biomass of individual trees to easily 
obtainable and nondestructive measurements, such as diameter, height and wood density (Chave 
et al., 2005). It has been demonstrated that choosing suitable allometric equations for each forest 
type is of great importance, because biomass and associated carbon estimates are highly sensitive 
to the choice of allometric equations (Jepsen, 2006). More accurate biomass estimations will be 
obtained when region-specific allometric equations will be developed (Pearson et al., 2007). For 
example, Iranmanesh et al. (2013) determined the allometric equations for different Persian oak 
trees (seed- and coppice-originated) by destructive method in Zagros forests, Iran. 

Zagros forests, with an area of 6x10° hm’, play important roles in water supply, soil 
conservation, climate change and socioeconomic balance in the semi-arid area, Iran 
(Sagheb-Talebi et al., 2014). However, Zagros forests significantly decreased due to forest 
clearance for cultivation, non-wood forest products exploitation, overharvesting and overgrazing, 
pollarding and firewood collection (Soltani et al., 2014). The Persian oak (Quercus brantii var. 
persica) is a dominant tree species of Zagros forests. We found three types of regeneration in oak 
based on their silvicultural characteristics, i.e., seed-originated oak, coppice oak and mixed 
(seed-originated and coppice) oak. Coppice stands that are characterized by sprout-clumps are 
formed asexually/vegetatively through sprouting. Seed-originated stands that account for less than 
10% in Zagros forests are dense and usually distribute in remote areas with old oak trees. And the 
diameter at breast height (DBH) of the old oak trees is higher than 80 cm. Mixed stands usually 
have two layers, in which the upper layer forms by seed-originated oak and the lower layer forms 
by sprout-clumps of oak (Sagheb-Talebi et al., 2014). Researches for assessing biomass and 
carbon stocks have been focused on Hyrcanian forests in Iran, however, less attention was paid to 
Zagros forests. This study analyzed the biomass and carbon stocks of three types of Persian oak to 
understand their contributions on regeneration of Zagros forests and provide a guideline for 
sustainable management of Zagros forests in Iran. 


2 Materials and methods 


2.1 Study area 


This study was conducted in Dalab valley with an area of 3000 hm? (33°41'01'-33°43/13''N, 
46°22'15"—46°25'27"E; 1200-2000 m a.s.l.) in Ilam Province, Iran. The study area has a cold and 
semi-arid climate with dry months extending from May to November. The mean annual 
precipitation is 525.0 mm (mainly occurs in winter) and the annual mean temperature is 16.9°C. 
In the study area, three types of Persian oak, i.e., seed-originated oak, coppice oak and mixed 
(seed-originated and coppice) oak were selected. There are no significant climatic differences in 
the study area among the three types of oak, and the topographic of slope, aspect, elevation and 
soil conditions are similar. Other accompanied tree species in the study area are Pistacia atlantica, 
Acer monspessulanum and Amygdalus eleagnifolia. 


2.2 Field measurements 


Thirty circular-shaped sample plots with an area of 1000 m? were assigned to the three types of 
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oak in a randomly systematic design. In each plot, indices of Persian oak including DBH (+5 cm), 
height, crown diameter and the number of sprouts (in coppice oak) were measured and recorded. 
All measured indices were calculated for each plot and averaged per hectare for each type of oak. 


2.3 AGB and AGCS 


For individual trees (DBH>5 cm) and coppice trees, AGB (kg/hm?) and AGCS (kg/hm?) were 
calculated using the Iranmanesh's allometric equations (Iranmanesh, 2013). 
Allometric equations for AGB and AGCS estimations of seed-originated oak were as follows: 
AGB=0.881X?-78, (1) 
AGCS=0.425X?:730, (2) 
where X is the average crown diameter (m). 
Allometric equations for AGB and AGCS estimations of coppice oak were as follows: 
AGB=2.534X?383, (3) 
AGCS=1,275X2262, (4) 
2.4 BGB and BGCS 


We estimated BGB of Persian oak based on a nondestructive method. We used the root to shoot 
ratio that was suggested by Askari et al. (2017) to estimate BGB. In this estimation, the average 
ratio of root to shoot in seed-originated and coppice oaks is considered as 0.8. Additionally, the 
coefficient of 0.5 for the conversion biomass into carbon was applied to obtain BGCS (Hernandez 
et al., 2004). 


2.5 LFB and LFCS 


To estimate the carbon content in tree litterfall, we randomly placed a quadrat with 1 mx1 m 
inside each plot (total 30 samples). Samples were dried at 70°C to a constant weight (Nunes et al., 
2013). 

The percentage of organic carbon (Ash) in each sample was measured by the combustion 
method in electric furnace. And carbon coefficient (C) of litterfall was estimated by the following 
equations (Allen et al., 1986): 


Ash = vH x 100% (5) 
W,- 1 f 

C=(100%-Ash)x0.58, (6) 

where Ash is the percentage of ash after burning (%), which represents the organic carbon; W3 is 

the sum of the weights of ash and vase (g); Wi is the weight of vase (g); W2 is the dry weight of 


sample (g); and C is the carbon coefficient (%). 
2.6 Statistical analysis 


The normality of the data was tested using Kolmogorov-Smirnov test. We used a one-way 
analysis of variance (ANOVA) to test the equality of means in AGB, BGB, AGCS and BGCS 
among the three types of oak. Then, a comparison of mean values of these variables was done by 
Duncan's test to find the differences. All statistical analyses were performed using IBM SPSS 
Statistics (version 16). 


3 Results 


3.1 Forest characteristics 


We selected 1036 trees among three types of oak and measured DBH, height, crown diameter and 
sprouts in 2017. We divided DBH and crown diameter into three diameter classes (i.e., <30, 
30-60 and >60 cm; <5, 5-15 and >15 m, respectively) and two height classes (i.e., <5 and 5-15 
m). Table 1 shows the quantitative characteristics of three types of oak. Results indicated that 
most of trees (93.70%) were in the smallest diameter classes (DBH<30 cm) and less than 1.00% 
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of trees had a DBH>60 cm (probably these trees had a seed origin) in coppice forest. Most of 
trees (61.03%) had a DBH of 30-60 cm and 33.10% of trees had a DBH<30 cm in 
seed-originated oak. Most of trees (65.48%) had a DBH<30 cm and 34.25% of trees had a DBH 
of 30-60 cm in mixed oak. From the above results, we found that coppice and mixed oaks had a 
higher percentage of trees with a DBH<30 cm, while seed-originated oak had the highest 
percentage of trees with a DBH>60 cm (Table 1). In coppice and seed-originated oaks, the highest 
frequency of tree height was observed in the class of 5-15 m, while in mixed oak, the highest 
frequency of tree height in the class of <5 m was observed (Table 1). 


Table 1 Diameter at breast height (DBH), height and crown diameter in the three types of Persian oak 


Type DBH(cm)  Gfwees  ofmees (06) TAMOD  ofuees  diameter(m) of trees 
Coppice oak <30 357 93.7 <5 161 <5 229 
30-60 21 5.5 5-15 220 5-15 152 
>60 3 0.8 - - >15 0 
Total 381 100.0 Total 381 Total 381 
Seed-originated oak <30 96 33.1 <5 12 <5 26 
30-60 177 61.03 5-15 278 5-15 261 
>60 17 5.87 - - >15 3 
Total 290 100.0 Total 290 Total 290 
Mixed oak <30 239 65.5 <5 200 <5 292 
30-60 125 34.2 5-15 165 5-15 73 
>60 1 0.3 - - >15 0 
Total 365 100.0 Total 365 Total 365 


Note: - means no value. 


3.2 Biomass and carbon stocks 


Table 2 shows the average biomass and carbon stocks in the three types of oak. The highest value 
in AGB was observed in seed-originated oak (76,043.25 kg/hm?) and followed by coppice oak 
(19,345.73 kg/hm?). In addition, the highest values of LFB and LFCS (3298.33 and 1520.48 
kg/hm?, respectively) were found in seed-originated oak. As shown in Table 2, the mean values of 
AGB (786.65 kg/hm?) and AGCS (380.05 kg/hm?) of single tree in seed-originated oak was 
higher than those of each group of sprout in coppice oak. Meanwhile, the lowest values of AGB 
(113.02 kg/hm°) and AGCS (53.71 kg/hm?) were found in mixed oak. 


Table 2 Biomass and carbon stocks in the three types of Persian oak 


Type T AGB AGCS BGB, BGCS LFB , LFCS 
(kg/hm?) (kg/hm?) (kg/hm?) (kg/hm?) (kg/hm?)  (kg/hm3 

Coppice oak Sprout 152.33? 73.81? 121.86" 60.93? - - 
Plot 1934.57? 937.32? 1547.65? 773.82 0.13? 0.06° 
Hectare 19,345.73? 9373.20? 15,476.58? 7738.29" 1334.00° 598.90° 

Seed-originated oak Tree 786.65? 380.05° 629.32? 314.66? e s 
Plot 7604.25" 3673.77 6083.46? 3041.70? 0.33? 0.10? 
Hectare 76,043.25? 36,737.79? 60,834.60° 30,417.30? 3298.33? 1520.48? 

Mixed oak Tree or sprouts 113.02 53.715 90.415 45.205 - - 
Plot 1375.00° 653.46° 1100.00° 550.00° 0.08° 0.04° 
Hectare 13,750.90° 6534.65° —-11,000.72° 5500.36° 812.33° 379.91° 


Note: AGB, aboveground biomass; AGCS, aboveground carbon stock; BGB, belowground biomass; BGCS, belowground carbon stock; 
LFB, litterfall biomass; LFCS, litterfall carbon stock. - means no value. Different lowercase letters indicate significant differences 
among the three types of oak at P<0.05 level. 


Figure 1 shows the percentages of AGCS in different DBH classes of the three types of oak. 
The results showed that trees with DBH<30 cm had the highest percentage of AGCS (81.45%) in 
coppice oak, while trees with 30-60 cm DBH had the highest percentage of AGCS in 
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Fig. 1 Percentage of AGCS (aboveground carbon stock) in different DBH (diameter at breast height) classes of 
the three types of oak. C, coppice oak; SO, seed-originated oak; M, mixed oak. 


seed-originated (71.1%) oak and followed by the mixed (63.2%) oak. Trees with DBH>60 cm 
contributed the least percentage of AGCS in the three types of oak. This result can be attributed to 
the lowest density of stems with the highest DBH class in the three types of oak. 

Based on the result of ANOVA, there were significant differences in DBH, height, crown 
diameter, LFCS, AGB, AGCS and BGCS among the three types of oak (Table 3). For all the 
variables, seed-originated oak had the significantly highest values than the other two types of oak 
(Fig. 2). The amounts of AGB, AGCS, BGCS and LFCS in coppice and mixed oaks were similar 
without significant differences (Fig. 2). 


Table 3 Results of ANOVA for the quantitative characteristics of oak forests 


Source of variation Sum of square df Mean of square F P 
Between groups 67,160.89 2 33,580.44 
DBH Within groups 154,547.61 1033 149.75 224.23 0.00™ 
Total 221,708.50 1035 
Between groups 2869.45 2 1434.72 
Height Within groups 4337.48 1033 4.19 341.69 0.00 
Total 7206.93 1035 
Between groups 3261.06 2 1630.53 
a, Within groups 4146.25 1033 4.01 406.23 0.00" 
Total 7407.32 1035 
Between groups 89,483,666.69 2 44,741,833.34 
AGB Within groups 177,789,569.70 1033 172,109.94 259.96 0.00 
Total 267,273,236.40 1035 
Between groups 20,956,004.11 2 10,478,002.05 
AGCS Within groups 41,771,163.30 1033 40,436.75 259.12 0.00 
Total 62,727,167.41 1035 
Between groups 10,854,014.11 2 5,427,007.05 
BGCS Within groups 31,761,162.30 1033 30,746.52 219.11 0.00 
Total 42,615,176.41 1035 
Between groups 20,000,198.00 2 10,990,000.00 
LFCS Within groups 8,101,614.98 1033 93,122,011.00 118.03 0.00 
Total 300,000,008.00 1035 


Note: DBH, diameter at breast height; AGB, aboveground biomass; AGCS, aboveground carbon stock; BGCS, belowground carbon 
stock; LFCS, litterfall carbon stock. ** indicates significant difference among the three types of oak at P<0.01 level. 


Figure 3 shows the amount of carbon stocks in different carbon pools of the three types of oak. 
In this study, carbon stocks included aboveground and belowground tree components and litterfall, 
i.e., AGCS, BGCS and LFCS. The maximum value of carbon stocks was found in seed-originated 
oak and followed by coppice oak. The minimum value of carbon stocks was found in mixed oak. 
In addition, AGB accounted for the highest percentage of carbon stocks, ranged from 52.64% in 


Ali MAHDAVI et al.: Biomass and carbon stocks in three types of Persian oak (Quercusi branitii... 


mixed oak to 53.49% in seed-originated oak. AGCS in coppice oak accounted for 52.93% of the 
total carbon stocks. In the belowground component, the percentages were almost equal for all the 
three types of oak. LFCS accounted for 2.21%, 3.06% and 3.38% of the total carbon stocks for 
seed-originated, mixed and coppice oaks, respectively. 
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Fig. 2 Mean values of variables of the three types of oak. DBH, diameter at breast height; AGB, aboveground 
biomass; AGCS, aboveground carbon stock; BGCS, belowground carbon stock; LFCS, litterfall carbon stock; C, 
coppice oak; SO, seed-originated oak; MF, mixed oak. Different lowercase letters indicate significant difference 
among the three types of oak at P<0.05 level. 
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Fig. 3 Carbon stocks in different carbon pools of the three types of oak. AGCS, aboveground carbon stock; 
BGCS, belowground carbon stock; LFCS, litterfall carbon stock. C, coppice oak; SO, seed-originated oak; MF, 
mixed oak. Different lowercase letters indicate significant difference among the three types of forest at P<0.05 
level. 


4 Discussion 


We observed three types (coppice, seed-originated and mixed oaks) of regeneration in Persian oak 
of Zagros forests in Dalab valley, Ilam Province, Iran. Therefore, knowledge of biomass and 
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carbon stocks of this species, especially in coppice oak (the dominant forest type) may play an 
important role in valuing the ecosystem and management planning of Zagros forests (Zianis et al., 
2005; Iranmanesh et al., 2014). Results indicate that biomass and carbon stocks were significant 
different among the three types of oak. In this study, the average biomass and carbon stocks of 
each tree or sprouts were 113.02 and 53.71 kg/hm? for mixed oak, 152.33 and 73.81 kg/hm? for 
coppice oak, and 786.65 and 380.05 kg/hm? for seed-originated oak, respectively, which had 
significant difference (P<0.05) among the three types of oak (Table 3). From the result, we can 
found that average biomass of trees in seed-originated oak was approximately five times higher 
than those of trees or sprouts in coppice and mixed oaks. Previous study found that AGB and 
AGCS were positively related to DBH, height and crown diameter of the tree (Ali et al., 2016). 
Therefore, this result is probably due to the presence of larger trees in seed-originated oak than 
sprouts in coppice oak (Table 1). The estimated average biomass and carbon stocks of spouts in 
coppice oak in this study were almost similar to the result of Iranmanesh et al. (2014) (146.3 and 
70.7 kg/hm?, respectively), but contrasted with the result of Alinejadi et al. (2016) (230.0 and 
115.0 kg/hm?, respectively). Moreover, the average biomass of single tree (786.65 kg/hm?) in 
seed-originated oak was higher than the estimations of Iranmanesh et al. (2014) and Alinejadi et 
al. (2016) (374.1 kg/hm? and 252.0 kg/hm/?, respectively). The reason for this discrepancy in 
biomass of seed-originated oak may be due to the variations in characteristics, structure, stand 
density, proportion of large trees, productivity and history of forests in these regions (Tran et al., 
2013; Mensah et al., 2016; Wang et al., 2016). 

Seed-originated oak had the highest values of AGB (76,043.25 kg/hm?) and AGCS (36,737.79 
kg/hm?), compared with other two stands (Fig. 3). Our findings are similar to the studies of 
Iranmanesh et al. (2014) and Alinejadi et al. (2016), who found that biomass and carbon stocks 
per hectare were higher in seed-originated oak than in coppice oak. Despite the similarity of 
results, the estimated biomass and carbon stocks of seed-originated oak in this study were 
approximately four times higher than those in Iranmanesh et al. (2014) and Alinejadi et al. (2016). 
The difference in AGB and AGCS among the three types of oak probably arise from forest 
structure, stand age and illegal wood cutting (Askari et al., 2017). Forest structure is original and 
identical in seed-originated stand in comparison with the other two types of oak. As shown in 
Table 1, the distribution of trees was almost uniform in the first two classes of DBH (i.e., <30 and 
30-60 cm) in seed-originated oak. Also crown diameter plays an important role in the total 
biomass of Persian oak trees and there are more trees with big crown diameter in seed-originated 
oak. However, the coppice shoots had different stems because of variations in their sprouting 
capacity. Crown diameter in the coppice shoots is smaller than that in seed-originated oak (Table 
1; Iranmanesh et al., 2014). Another important factor for biomass partitioning within a tree 
species is tree age. With increasing age, the proportion of stem in the biomass becomes more 
important. Then proportion of canopy biomass increases with increasing tree size. Thus, forest 
intervention and transformation from original forest to destroyed forest are the main factors that 
cause changes in forest biomass (Peichl and Arain, 2007; Sanquetta et al., 2011). 

The results of the study showed that litterfall contained carbon stocks of 379.91—1520.48 
kg/hm? (Table 2; Fig. 3). Litter layer is a significant reservoir of carbon in forest ecosystems 
especially in soil organic carbon stock (Lorenz and Lal, 2010). Litter biomass and carbon stocks 
in seed-originated oak is significantly higher than other two types of oak because of the more 
number of trees per hectare and the wider crown of trees. Litter layer, which result from the 
annual amount of litterfall accumulation, plays an important role in forest carbon budgets (Lorenz 
and Lal, 2010). According to previous study, forest litter and soils store approximately 3-fold 
more carbon than living plants (Wang et al., 1999), and carbon pool of litter layer may store 
5%-20% of the total stand carbon storage (Karjalainen, 1996). 

Our findings showed that total carbon sequestration (68,675.57 kg/hm?) in seed-originated oak 
was significantly higher than the other two types of oak. A similar trend to the reports of 
Iranmanesh et al. (2014) and Askari et al. (2017) was found in Zagros forests. The highly variable 
of the total biomass and carbon stocks in studied forests is supported by the strong role of 
structural variable of forests, which could also be used to management practices and conservation 
plans. 
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5 Conclusions 


In this study, we estimated the biomass and carbon stocks of aboveground and belowground tree 
components and litterfall among the three types of Persian oak. As expected, the origins of trees 
(seed or coppice) and forest structural variables including DBH, stand density and proportion of 
large trees in forests were found to explain a high variability of the estimated biomass and carbon 
stocks in studied forest types. The differences in biomass and carbon stocks in Zagros forests 
indicate the need for site-specific biomass and carbon stocks equations for the accurate 
assessment in these areas. Our result can help the managers and decision-makers to consider the 
role of Zagros forests in reducing atmospheric carbon dioxide and its carbon storage potential. 
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